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^ . Abstract 

From about 4 million hadronic Z° decays recorded by the OPAL detector on and near 
to the Z° resonance, we select a sample of more than 570 000 inclusively reconstructed 
B mesons. Orbitally-excited mesons Bj are reconstructed using TStt^ combinations. In- 
dependently, B* mesons are reconstructed using the decay B* — * B7. The selected B* 
candidates are used to obtain samples enriched or depleted in the decay B} —* B*7r :l= (A), 
where (X) refers to decay modes with or without additional accompanying decay parti- 
cles. From the number of signal candidates in the mass spectra of these two samples, 
we perform the first measurement of the branching ratio of orbitally-excited B mesons 
decaying into B*7r(X): 

X; BR(B} -» B*7TpQ) = 0.85 loff ±0.12, 

where the first error is statistical and the second systematic. If B j decay modes other than 
single pion transitions can be neglected the measured ratio corresponds to the branching 
ratio BR(B} -> B*vr). 

In the framework of Heavy Quark Symmetry, a simultaneous fit to the B7r ± mass 
spectra of the samples enriched or depleted in Bj — > B*7r =l= (X) decays yields the mass and 
the width of the Bi(3/2) state, as well as the branching ratio of B} mesons decaying into 
B*tt: 

M(Bi(3/2)) = ( 5.738 ± g;ggg ±0.007) GeV/c 2 
r(B!(3/2)) = (18lJt;i)MeV/c 2 
BR(B}-B*tt) = 0J4±° it±° il , 

where the uncertainties are statistical and systematic, respectively. 



(Submitted to Eur. Phys. J.) 
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1 Introduction 

An important prediction of Heavy Quark Effective Theory (HQET) is the existence of an ap- 
proximate spin-flavour symmetry for hadrons containing one heavy quark Q (mg ^> Aqcd) 0- 
In the limit mg — > oo, mesons composed of a heavy quark Q and a light quark q are charac- 
terised by the spin of the heavy quark, Sq, the total angular momentum of the light quark, 
jq = S q + L, and the total angular momentum, J, where S q and L denote spin and orbital an- 
gular momentum of the light quark, respectively. In the heavy quark limit, both Sq and j q are 
good quantum numbers and the total angular momentum of the meson is given by J = Sq + j q . 
For L — 1, there are four states with spin-parity J p = + , 1 + , 1 + and 2 + . If the heavy quark Q 
is a bottom quark, these states are labelled Bq, Bi for both 1 + states[] and B 2 0, respectively. 
The four states, commonly called B}, or alternatively B** 0, are grouped into two sets of de- 
generate doublets, corresponding to = 1/2 and j q = 3/2 as indicated in Table |Tj. Parity and 
angular momentum conservation put restrictions on the strong decays of these states to B^*'tt f\ 
(see Figure |TJ). The + state can only decay to Btt via an S-wave transition, the lw 2 to B*7r 
via an S-wave transition, the ljj/ 2 to B*7r via a D-wave transition, and the 2 + state can decay 
to both Btt and B*tt via D-wave transitions only. States decaying via an S-wave transition are 
expected to be much broader than the states decaying via a D-wave transition ||. In addition 
to the single pion transitions, decays to B*tttt and Btttt are also possible. In the case of di-pion 
transitions, all four B} states are allowed to decay to B* as well as to B. Although these decays 
are phase-space suppressed, intermediate states with large width like B} — > B^*'p — > B^tttt 
may cause a significant enhancement of the B*titt and Btttt final states ||. Additional B} decay 
modes with other than one or two accompanying pions are expected to be strongly suppressed 
but can not be excluded. Therefore, the notations B*7r(X) and B7r(X) are chosen to refer to 
the final states of B} decays. 

Given the HQET predictions listed in Table [|, the four B} states are expected to overlap 
in mass. So far, in analyses from LEP experiments and from CDF || B} mesons are 

reconstructed in the Btt final state only, observing one single peak in the Btt mass spectrum. 
This is not sufficient to resolve any substructure of the four expected B} states. In addition, 
for decays to B*tt where the photon in the decay B* — > B7 is not detected, the reconstructed 



Btt mass is shifted by Mb — M-q* = —46 MeV/c 2 . A recent analysis [ 10 1 tries to cope with these 
problems by constraining all properties of the four B} states according to HQET predictions 
except for the masses and widths of Bi(l/2) and B 2 . 

In this paper a different approach is presented. Using information from the photon in the de- 
cay B* — > B7, the B} — > B*7r ± (A) decays are statistically separated from the B} — > B7r ± (A) de- 
cays. This allows a model-independent measurement of the branching ratio BR( B} — ► B*tt(X)). 
Assuming that B^7rX decays produce small contributions to the B} width, this method gives 
insight into the decomposition of the B} into the states allowed to decay to B-zr (Bq and B 2 ) 
from the other states that can only decay to B*tt. 

The Bit invariant mass spectrum is also fit in the context of HQET expectations. This 
requires high statistics, extensive background studies and the tagging of the B* decay. Due to 
the complexity of the B} signal and its different decay modes, several constraints provided by 



1 In the limit toq — > 00, the notations Bi(l/2) and Bi(3/2) are used. In the case of mixing of the J = 1 
states, the notations Bi(H) and Bi(L) arc used to distinguish the physical states. 

2 Throughout this paper, we use the Particle Data Group notation Bj for orbitally-excited B mesons. 

3 Throughout this paper, B^*^7r denotes the final states Btt and B*tt. The notations B^tttt and B*tt(X) are 
to be interpreted in the same way. 
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\r* \t 1 21 
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°l/2 


5.738 
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Bi 


1 + 
-"-3/2 
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0.020 


B* 


2+ 

^3/2 


5.733 


0.025 


(B*7r)D-wave, 


5.737 


0.025 








(B7r) D _ wave 







Table 1: Masses, widths and dominant decay modes based on theoretical predictions P,pl|,|TT|- 
D^] and the corresponding Monte Carlo input values used in the analysis. Recent calculations 
using a bag model predict widths of T(B*) = 0.141 GeV/c 2 and r(Bi(l/2)) = 0.139 GeV/c 2 for 
the broad states 



HQET are imposed. From the fitting procedure, we obtain model-dependent measurements of 
the mass and width of B!(3/2) as well as BR(B} — > B*7r). 

The paper is organised as follows: the next section describes the data sample and the event 
simulation. In Section [| the analysis method is presented. The B reconstruction is described 
in Section Section [5] contains the photon reconstruction. The pion reconstruction and the 
total Bit mass spectrum are presented in Section |6|. The BR(B} — ► B*7i(X)) measurement and 
results from a simultaneous fit to the Bti^ 1 mass spectra of the samples enriched or depleted in 
the decay B} — > B*tc ± (X) are presented in Section [7|. Systematic uncertainties are evaluated 
in Section || A discussion of the results and conclusions are given in Section ^. 



2 Data sample and event simulation 

The data used for this analysis were collected from e + e~ collisions at LEP during 1991-1995, 
with centre-of-mass energies at and around the peak of the Z° resonance. The data correspond 
to an integrated luminosity of about 140 pb _1 . A detailed description of the OPAL detector 
can be found elsewhere |15],|16|. 



Hadronic events are selected as described in [17] > giving a hadronic Z selection efficiency of 



(98.4 ± 0.4) % and a background of less than 0.2 %. A data sample of about 4 million hadronic 
events is selected. Each event is divided into two hemispheres by the plane perpendicular to 
the thrust axis and containing the interaction point of the event. The thrust axis is calcu- 
lated using tracks and electromagnetic clusters not associated with any tracks. To select events 
within the fiducial acceptance of the silicon microvertex detector and the barrel electromag- 
netic calorimeter, the thrust axis direction?] is required to satisfy | cos 6^1 < 0.8. Monte Carlo 
simulated samples of inclusive hadronic Z° decays are used to evaluate efficiencies and back- 
grounds. The JETSET 7.4 parton shower Monte Carlo generator ||18||, with parameters tuned 



by OPAL |I9| and with the fragmentation function of Peterson et al. |20[ for heavy quarks is 
used to generate samples of approximately 10 million hadronic Z° decays, 2 million Z° — > cc and 
5 million Z° — > bb decays. The generated events are passed through a program that simulates 



the response of the OPAL detector |21[ before applying the same reconstruction algorithms as 

4 In the OPAL coordinate system, the x axis points towards the centre of the LEP ring, the y axis points 
upwards and the z axis points in the direction of the electron beam. 9 and cf> are the polar and azimuthal angles, 
and the origin is taken to be the centre of the detector. 
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for data. All generated Monte Carlo samples contain L = 1 states for bottom and charmed 
mesons, as well as vector meson partners of the ground states. The generated production rates, 
masses and widths of all resonant states are consistent with experimental measurements when 
available and with theoretical predictions elsewhere (see also Table |l|). 

3 Analysis overview 

The analysis is based on the reconstruction of B* in the B7 final state and a separate recon- 
struction of B} in the B7r ± final state. A direct reconstruction of B} decaying to B*7r, B* — > B7 
giving B77r ± in the final state is inappropriate because of the large combinatorial background 
and the insufficient detector resolution. Therefore, our approach employs a statistical separa- 
tion of B} -> B*vr ± (X) from B} -> Bvr ± (X) decays. 

B mesons produced in Z° — > bb events are reconstructed inclusively to achieve high efficiency. 
No attempt is made to reconstruct specific B decay channels. On the contrary, properties 
common to all weakly decaying b hadrons are used for the B reconstruction. For each B 
candidate, a weight W(B*) is formed where W(B*) represents the probability for the B to have 
come from a B*. The probability W(B*) is based on the reconstruction of photon conversions 
and of photons detected in the electromagnetic calorimeter. All B candidates are then combined 
with charged pions to form B} meson candidates. Using the weight W(B*), we derive two 
mutually exclusive subsamples of B^ combinations, one enriched and the other depleted in 
its B* content. Invariant Btt mass distributions are formed for both samples. The shape 
of the non-B} background of the two distributions is taken from Monte Carlo simulation and 
normalised to the data in the upper sideband region and subtracted from the corresponding 
data distributions. The branching ratio BR(B} — > B*7r(X)) is obtained from the observed 
number of B} and the different efficiencies for B} — > B*tt ± (X) and B*j — > B7r ± (X) decays in 
the B*-enriched and the B*-depleted samples. Applying a simultaneous fit to the Bn mass 
spectra of both samples several details of the B} four-state composition and of the B} decay 
modes are extracted. 

Whereas the BR(B} — ► B*7r(X)) result obtained from counting the number of B} signal 
entries of the samples enriched or depleted in the decay B} — > B*7r ± (A) is model-independent 
and does not rely on the shape of the Bn signal, the fit to the Bit mass spectra makes use of 
HQET assumptions on the composition and the decay modes of the B} signal. 



4 Selection and reconstruction of B mesons 

B mesons are reconstructed using an extended version of the method used in earlier analyses 



2!|. Since the reconstructed B mesons are used to form B* and B} candidates, the B 



reconstruction is tuned to minimise the uncertainties on the B direction and energy, while 
maintaining a high reconstruction efficiency. 

4.1 Tagging of Z° — ► bb events 

To achieve optimal b-tagging performance, each event is forced to a 2-jet topology using the 
Durham jet-finding scheme ||23|| . In calculating the visible energies and momenta of the event 
and of individual jets, corrections are applied to prevent double counting of energy in the case 
of tracks with associated clusters [24]]. A b-tagging algorithm is applied to each jet using three 
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independent methods: lifetime tag, high pr lepton tag and jet-shape tag. A detailed description 
of the algorithm can be found in . The b-tagging discriminants calculated for each of the 
jets in the event are combined to yield an event b likelihood B event . For each event, B event > 0.6 
is required. After this cut, the Z° — > bb event purity is about 96%. The cut on the direction 
of the event thrust axis, | cos^tI < 0.8, as described in Section 0, removes roughly a quarter of 
all Z° — > bb events and after the cut on B cvent , the total b event tagging efficiency with respect 
to all produced Z° — > bb events is about 49%, where these numbers are obtained from Monte 
Carlo simulation. At this stage, about 750 000 b hadron candidates are selected. 



4.2 Reconstruction of B energy and direction 

The primary event vertex is reconstructed using the tracks in the event constrained to the av- 
erage position and effective spread of the e + e~ collision point. For the b hadron reconstruction, 
tracks and electromagnetic calorimeter clusters with no associated track are combined into jets 



using a cone algorithm [£6[ with a cone half-angle of 0.65 rad and a minimum jet energy of 
5.0 GeV (]. The two most energetic jets of each event are assumed to contain the b hadrons. 
For each jet we reconstruct the energy and direction. 

In each hemisphere defined by the jet axis, a weight is assigned to each track and each cluster, 
where the weight corresponds to the probability that any one track or cluster is a product of 
the b hadron decay. The b hadron is reconstructed by summing the weighted momenta of the 
tracks and clusters. The reconstruction algorithm is applied to all b hadron species and is 100% 
efficient. Since this analysis aims at the reconstruction of _B Uj d mesons which make up about 
80% of the b hadron sample, b hadron candidates are referred to as B mesons in the following. 
Details of the reconstruction method are provided below. 



4.2.1 Calculation of track weights 

Two different types of weights are assigned to each track: 

• coVtx, calculated from the impact parameter significances of the track with respect to both 
the primary and secondary vertices; 

• co>nn, the output of a neural network based on kinematics and track impact parameters 
with respect to the primary vertex. 

The calculation of u vtx requires the existence of a secondary vertex, whereas u;nn does not and 
is therefore available for all tracks. The search for detached secondary vertices proceeds as 
follows: 

Each jet is searched for secondary vertices using a vertexing algorithm similar to that de- 
scribed in m, making use of the tracking information in both the r — <p and r — z planes if 
available. If a secondary vertex is found, the primary vertex is re-fitted excluding the tracks 
assigned to the secondary vertex. Secondary vertex candidates are accepted and called 'good' 
secondary vertices if they contain at least three tracks and have a decay length greater than 
0.2 mm. If there is more than one good secondary vertex attached to a jet, the vertex with the 
largest number of significant^ tracks is taken. If there is a tie, the secondary vertex with the 

5 The cone jet-finder provides the best b hadron energy and direction resolution compared to other jet finders 
studied here. 

A track is called significant if its impact parameter significance with respect to the primary vertex is larger 
than 2.5. 
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larger separation significance with respect to the primary vertex is taken. If a good secondary 
vertex is determined, a weight is calculated for each track in the hemisphere of the jet using 
the impact parameter significance of the track with respect to both the primary and secondary 
vertices. This weight is given by 

R(Wv) m 

vtx R(b/r,)+ R(d/a) ' [ ) 

where b and r] are the impact parameter and its error with respect to the secondary vertex, and 
d and a are the same quantities with respect to the primary vertex. R is a symmetric function 
describing the impact parameter significance distribution with respect to a fitted vertex. The 
u; v tx distribution for tracks of hemispheres with a good secondary vertex is shown in Figure [|a 
and compared with the corresponding Monte Carlo distribution. The weight u vtx shows a weak 
correlation with the momentum of the track. 

For each track, the weight wnn is calculated using an artificial neural network trained to 
discriminate b hadron decay products from fragmentation tracks in a jet. The neural network 
was trained using as inputs the scaled track momentum x p = p/E heam , the track rapidity 
relative to the estimated B direction, the impact parameters of the track with respect to the 
primary vertex in the r — and r — z planes and the corresponding errors on the impact 



parameters [p8 |. As a preliminary estimate, the jet axis is taken as the estimated B direction. 
The unn distribution is shown in Figure |||b. If a good secondary vertex exists, the track weight 
cunn is combined with the vertex weight u vtx using the prescription 

^NN ' ^vtx /ri N 
(1 — Wnn) • (1 — ^vtx) + ^NN • ^vtx 

The weight Ut r in Equation is approximately the probability that the track is a b hadron 
decay product. In the case where there is no good secondary vertex in the jet, the total track 
weight u tT is simply given by u tT = c^nn- The combined weight u tT for tracks of all hemispheres 
is shown in Figure 0c. 



4.2.2 Calculation of cluster weights 

Similar weights are calculated for energy clusters reconstructed in the electromagnetic and 
hadronic calorimeters to represent the probability the clusters came from a B hadron decay. 
Weights u) ec \ and uvi are assigned to each electromagnetic and hadronic cluster in the hemi- 
sphere of the B meson based on their rapidity with respect to the estimated B direction. The 
weight is equal to the probability, calculated using the Monte Carlo simulation as a function 
of the cluster energy, that the cluster came from the decay of a B meson. Clusters associated 
with a track have the estimated energy of the track subtracted. 



4.2.3 Calculation of B direction 

The B momentum is calculated iteratively by a weighted sum of all tracks and clusters in the 
hemisphere: 

-^track -^ccal -^hcal 

V = ' Pi + U ecL,i -Pi+ ^hcl,i " Pi (3) 

i=l i=l i=l 

where A^ track , iV eca i and -/V hcal denote the number of tracks, electromagnetic clusters and hadronic 
clusters, respectively. The rapidity calculation, for both tracks and clusters, is initially per- 
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formed relative to an estimate of the B meson direction]]. The weights are then recalculated 
with the rapidity determined using the new B direction estimate. 

An estimate of the B hadron direction is made for all B candidates based on the weighted 
momentum sum of the tracks and clusters in the jet. In addition, the vector from the primary 
vertex to the secondary vertex yields an estimate of the B direction for those jets where a good 
secondary vertex has been identified. When both estimates are available the weighted average 
is taken, using the calculated uncertainties of each direction estimate. The covariance matrices 
of the primary and secondary vertices determine the error on the B flight direction. The error 
on the momentum sum is estimated by removing each term in turn from the sum in Equation |^, 
calculating the change in the B direction caused by this omission and adding up in quadrature 
the corresponding error contributions from each track and cluster. The final estimate of the B 
direction is obtained by taking the error-weighted sum of the B direction calculated with the 
momentum sum method and the B direction obtained from the primary and secondary vertex 
positions. The direction information in the r — z plane of the secondary vertex is only used if 
the vertex is built with tracks that left at least four hits in the ^-layers of the silicon microvertex 
detector (the maximum number of these hits per track is two). 

The error Aa on the weighted sum of both B direction estimators described in the previous 
paragraph is a measure for the quality of the B direction^. To improve the resolution on the 
B direction, which in turn dominates the Bit mass resolution, a cut on Aa is imposed. Since 
this analysis aims at a separation of some of the B} states by reconstructing different B} decay 
channels rather than obtaining a very good Bit mass resolution, the cut Aa < 0.035 is rather 
loose. This cut removes the 20% of the B candidates with the poorest direction resolution, 
mainly those with no associated good secondary vertex. 



4.2.4 Calculation of B energy 

The resolution on the total energy of the B candidate can be significantly improved by con- 
straining the total centre-of-mass energy, Ecu, to twice the LEP beam energy. Assuming a 
two-body decay of the Z°, we obtain: 

p E CM ~ ^recoil + M B ( A\ 

B = ^ > ( 4 ) 

where the mass of the b hadron is set to the B meson mass Mb = 5.279 GeV/c 2 and M reco ii 
denotes the mass recoiling against the B meson. The recoil mass and the recoil energy -Erccoii 
are calculated by summing over all tracks and clusters^ of the event weighted by (1 — Uj) and 
assuming the particle masses used in the calculation of Ei. To account for the amount of 
undetected energy mainly due to the presence of neutrinos, the recoil mass is scaled by the 
ratio of the expected energy in the recoil to the energy actually measured: 

M - M EcM ~ Eb (K\ 

-'"recoil, new -'"recoil, old ' ,-, \i'J 

-^recoil 

where E& is taken from Equation [|. The new recoil mass value M recoi i new obtained from 
Equation || is substituted into Equation f| and the calculation of E B is iterated. After two 

The initial input for this axis is the jet direction calculated using tracks and unassociated electromagnetic 
clusters. 

8 In the case where no good secondary vertex exists, Aa is simply given by the uncertainty on the momentum 
sum. 

9 Tracks and clusters not contained in the hemisphere of the B meson candidate have weights u>i = 0. u>i 
denotes the weight Wtiyt, w cci,i and u>hci,i for tracks, electromagnetic clusters and hadronic clusters, respectively. 
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iterations the uncertainty on the B meson energy is minimised. A minimum B energy of 
15 GeV is required to further improve the energy resolution. 

After all these cuts, the narrower Gaussian from a two Gaussian fit to the difference between 
the reconstructed and generated B meson energy has a = 2.3 GeV, and 86% of the entries are 
contained within 3a. The distribution of the difference between the reconstructed and generated 
cf) angle of simulated B mesons can be described by a similar fit. The standard deviation of 
the narrower Gaussian is 14.2 mrad and 88% of the entries lie within 3a. The corresponding 
quantities describing the 9 resolution are a = 15.0 mrad and 89%, respectively. 

The complete B meson selection applied to the full data sample results in 574 288 tagged 
jets with a b purity of about 96%, as estimated from Monte Carlo. About 75% of the selected 
jets contain a good secondary vertex. 



5 The decay B* — ► B7 

The photon produced in the decay B* — > B7 has an energy of about 46 MeV in the rest frame 
of the B*. The mean energy of the photon in the laboratory frame is approximately 350 MeV, 
with a maximum energy below 800 MeV. Due to the kinematics of the process, these photons 
are produced predominantly in the core of the jet. The high particle density in this region gives 
rise to a high background level when identifying the photon. Since a high B* reconstruction 
efficiency is crucial for this analysis, photons are reconstructed in two ways: from energy 
deposits in the electromagnetic calorimeter and from converted photons in the tracking volume. 
The conversion probability within the OPAL tracking system for photons coming from the decay 
B* — > B7 is approximately 8%. 



5.1 Reconstruction of photon conversions 

The reconstruction of converted photons is optimised for the low energy region. The selection 
algorithm is partially based on quantities that have been used in earlier analyses |29| but tuned 



to obtain high efficiency rather than very good angular and momentum resolution. Given the 
low energy carried by these photons, we ignore calorimetry information and only use tracking 
information for the reconstruction of converted photons. 

Tracks with a total momentum p below 1.0 GeV/c, opposite charge assignment and a mea- 
sured dE/dx within three standard deviations of the expected value for electrons are combined 
into pairs. For each pair, the track with the greater scalar momentum is required to have a 
transverse momentum p t > 50 MeV/c with respect to the beam axis and at least 20 hits in 
the central jet chamber. For the track with lower momentum, a minimum p t of 20 MeV/c is 
required. The asymmetric selection cuts for the two tracks in a pair guarantee at least one well 
measured track and reflect the fact that the electron and the positron of a converted photon 
tend to have different momenta in the laboratory frame. To suppress random track combina- 
tions, the distance of closest approach between the two tracks of a pair in the r — <p plane has 
to be smaller than 1.0 cm with an opening angle between the tracks at their point of closest 
approach smaller than 1.0 rad. 

In order to make optimal use of all the available information, the following measured quan- 
tities for each conversion track pair candidate are fed into a neural network: 

• the distance of closest approach between the two tracks in the r — (f> plane; 
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• the radial distance with respect to the z axis of the first and last measured hits in the 
inner tracking chambers for each track; 

• the radial distance with respect to the z axis of the common vertex^ of both tracks 
obtained from a fit in the r — plane; 

• the impact parameter with respect to the primary vertex in the r — <fi plane of the recon- 
structed photon; 

• the invariant mass of the track pair assuming both tracks to be electrons; 

• the transverse momentum relative to the z axis of the lower momentum track. 

All conversion candidates with a neural network output greater than 0.7 and a photon energy 
below 1.5 GeV are called 'good' conversion candidates for a given B meson candidate if the 
opening angle between the reconstructed B momentum vector and the reconstructed photon 
momentum vector is smaller than 90°. At this stage, an average of 0.82 good conversion 
candidates is selected per B candidate in both data and Monte Carlo. The candidate multiplicity 
distributions are shown in Figure 0a. The total efficiency to detect photons from the decay 
B* — ► B7 with the conversion algorithm is estimated from simulation to be (2.70 ± 0.01 s t a t)%- 
The efficiency is rather independent of the photon energy from 1.0 GeV down to 200 MeV where 
it rapidly drops to zero due to track selection requirements. The amount of fake conversions in 
the selected sample is estimated from Monte Carlo simulation to be (11.75 ± 0.04 stat )%. 

Fits to the difference between the reconstructed and generated photon energy in Monte 
Carlo are made using the sum of two Gaussians, both constrained to the same mean value. 
The narrower Gaussian has a standard deviation of 5 MeV at an energy of 200 MeV, and rises 
to 13 MeV at an energy of 750 MeV, and about 70% of the entries are contained within 3cr. 
Similar fits to the <fi and 9 resolutions give values of 3.4 mrad (70%) and 5.4 mrad (61%), 
respectively. 



5.2 Reconstruction of photons in the electromagnetic calorimeter 

Photons are also detected as showers in the barrel region of the electromagnetic calorimeter. 
The location and energy of these showers are obtained from a fit to the energy deposits in the 
individual lead glass blocks not associated with any track. The whole reconstruction method 
has been shown to work in the dense environment of hadronic jets down to photon energies as 
low as 150 MeV. The details of the reconstruction are given in ||30|| . 

Showers in the electromagnetic calorimeter are accepted as photon candidates if they have 
an energy in the range 200 MeV to 850 MeV and a photon probability P 7 > 0.20, where P 7 



is the output of a simplified neural network [30]. If the opening angle between such a shower 



and a reconstructed B candidate is less than 90°, this shower is considered a 'good' photon 
candidate for the corresponding B candidate. On average, there are 4.59 (4.38) good calorimeter 
photon candidates per B candidate selected in the data (Monte Carlo) sample. To correct for 
the observed discrepancy, the Monte Carlo is reweighted to the data distribution shown in 
Figure |3|b. The efficiency to detect a photon from the decay B* — > B7 in the electromagnetic 
calorimeter is estimated to be (14.52 ± 0.03 stat )% using Monte Carlo simulated events. The 
fraction of fake photons arising from tracks and neutral hadrons in the sample ranges from 



10 The z position of this vertex is fitted independently and the reconstructed photon vector is constrained to 
the z coordinate of the primary vertex to improve the accuracy of the 6 determination. 
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32% at a photon energy of 850 MeV up to 43% at a photon energy of 200 MeV. If compared 
with the selected conversion sample, the selection of B* — > B7 photons in the electromagnetic 
calorimeter has a much higher efficiency but lower purity. 

As with the converted photons, the energy resolution has been determined from Monte Carlo 
simulation using a double Gaussian fit. The narrower Gaussian has a width of 20 MeV at a 
photon energy of 250 MeV and increases up to 86 MeV at an energy of 800 MeV, and about 75% 
of the entries are contained within 3cr. Similar fits to the <ft and 9 resolutions give values of 3.6 
mrad (65%) and 3.6 mrad (72%), respectively. In contrast to the conversion sample, photons 
reconstructed in the electromagnetic calorimeter have much higher energy uncertainties, but a 
better 9 resolution. 

5.3 Reconstruction of B* — > B7 decays 

Each reconstructed B meson candidate is combined with all good conversion and calorimeter 
candidates to reconstruct B* candidates. The invariant mass of a B7 combination is defined as 



where Mb is 5.279 GeV/c 2 and a is the measured angle between the B meson and the photon 
candidate. The mass difference AM = Mb 7 — Mb between the B* candidate and the B is 
calculated by simply subtracting the nominal B mass of Mb = 5.279 GeV/c 2 from Mb 7 . 

The mass difference distributions of the conversion sample observed in the data and the 
corresponding Monte Carlo background are shown in Figure |]a. The background is normalised 
to the data in the sideband region 0.09 GeV/c 2 < AM < 0.20 GeV/c 2 . The background 
subtracted signal of Figure £|b is fitted to the sum of two Gaussians fixed to the same mean, 
where one of the Gaussians is allowed to have asymmetric width. The observed asymmetry of 
the mass resolution of the conversion sample is well simulated in the Monte Carlo and is due 
to the very loose track requirements of the lower momentum track of the conversion pair. A 
mass difference of AM = (45.87 ± 0.25 sta t) MeV/c 2 is obtained from the fit to the data, where 
the error is statistical only. This result agrees well with the current world average value of 
(45.78 ±0.35) Me V/c 2 @. 

The AM distribution using calorimeter photons is shown in Figure [|a. The background 
is taken from Monte Carlo simulation and normalised to the data in the sideband region 
0.10 GeV/c 2 < AM < 0.20 GeV/c 2 . The same fit function as for the photon conversion sample 
is used to obtain the mass difference AM from the background subtracted signal distribution 
in Figure |5|b. A value of (47.30 ± 0.61 sta t) MeV/c 2 is obtained from this fit, which is consistent 
within two sigma with the result from the fit to the conversion signal. 

For the B* sample reconstructed with photon conversions, the mass resolution is dominated 
by the uncertainty on the reconstructed B direction. For the calorimeter photon sample, the 
B* mass resolution suffers in addition from the energy resolution of the calorimeter. Due 
to the high background of fake photons and the moderate energy resolution at low photon 
energies, the signal-to-background ratio is rather poor for calorimeter photons. Therefore, 
uncertainties in the B* reconstruction using these photons are dominated by systematic errors 
on the background shape and energy calibration. All systematic uncertainties arising from the 
B* reconstruction will be discussed in Section ||. 




(6) 
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5.4 The B* probability W(B*) 

To select samples enhanced and depleted in B* mesons, a B* probability is assigned to each B 
candidate. This probability combines information from both conversion and calorimeter photon 
candidates and represents the probability that a B candidate is the true daughter of a B* meson. 
Only the best conversion and best calorimeter candidate assigned to any one B candidate are 
considered in the calculation of this probability, where the best candidate is defined as that 
which gives AM = M By — M B closest to the world average of 45.78 MeV/c 2 |§. 

This weight is constructed by parametrising the purity of the mass difference distribution in 
several variables in Monte Carlo simulation. For calorimeter photon candidates, this parametri- 
sation is performed as a function of the photon probability, P 7 (see Section ^2), and the total 



number of good calorimeter photon candidates found per B candidate. For each B candidate, 
a single weight is calculated by taking the simple mean of the weight resulting from each of the 
above parametrisations. 

Similarly, for conversion photon candidates, the parametrisation is performed in AM as 
a function of the total number of conversion candidates, and a weight is extracted as for the 
calorimeter candidates. The two weights obtained from conversion and calorimeter photons are 
combined by taking their mean. 

The resulting weight W(B*) is shown in Figure |6|a for Monte Carlo and data, and the 
contributions from jets containing a B* and jets containing no B* as seen in the simulation are 
shown. The primary features of the W(B*) distribution are: 

• a peak at W(B*) = 0.625, corresponding to B candidates with no associated good con- 
version or calorimeter photon candidate; 

• a peak at W(B*) = 0.632, containing B candidates with no good conversion candidate 
and a best calorimeter candidate having a B7 mass far away from the nominal B* mass; 

• a peak around W(B*) = 0.656, containing B candidates with the best calorimeter candi- 
date close to the nominal B* mass; 

• a peak at W(B*) = 0.715, containing B candidates with the best calorimeter candidate 
being close to the nominal B* mass and having a high photon probability P 7 ; 

• the tail towards high B* probabilities is made up by best conversion candidates very close 
to the nominal B* mass. 

The assignment of the specific photon candidate samples to the peaks and to the tail of the 
W(B*) distribution is based on Monte Carlo information. Figure |6|b shows the ratio e(B*)/e(B) 
versus W(B*). £:(B*) refers to the efficiency to select a B meson from a true B* — > B7 decay 
and e(B) is the efficiency to select a B meson which has not come from a B*. In general, the 
Monte Carlo simulation of W(B*) describes the data adequately. A comparison of the Monte 
Carlo and data distributions yields a x 2 P er degree of freedom of two. A cut on W(B*) allows 
the production of samples of B candidates with different B* fractions. Further details and 
systematic studies concerning W(B*) are given in Section R. 



6 Reconstruction of orbitally-excited B mesons 

All B} candidates, even those expected to decay into B*^, are reconstructed using the mea- 
sured four-momenta of the B meson and the pion. B candidates are selected and reconstructed 
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as described in Section [| and combined with charged pion candidates. Pions produced in the 
decay of a B} will be referred to as 'signal pions'. Since the B} decays strongly, signal pions 
are expected to be associated to the primary event vertex rather than to a possible secondary 
vertex. In comparison with other pions created in the fragmentation process, signal pions are 
expected to have a large longitudinal momentum p\ with respect to the jet axis. These are the 
basic characteristics used to separate signal pions from B decay products and from fragmen- 
tation tracks. A significant number of non-resonant fragmentation pions are expected to be 
produced near a B meson. The kinematics of these pions is similar to the signal pions, giving 
rise to a combinatorial background in the invariant mass of Btt candidates. Tracks from B 
decay also contribute to the background due to the inability to unambiguously associate all B 
decay tracks with the secondary vertex. 



6.1 Pion selection 

The signal pion selection for this analysis makes use of techniques used in || and || . All tracks 
that are well measured according to a standard track selection |HJ are considered as possible 
signal pion candidates if they belong to the same jet as the B candidate. Additionally, the 
following selection cuts are applied in the given order: 

• The measured ionisation energy loss dE/dx has to be consistent with the expected value 
for pions within 2.6 standard deviations, if dE/dx information is available for this track. 



To suppress B decay tracks, the track weight wnn as described in Section fO| has to be 
smaller than 0.9. 

The B decay track rejection is improved by the requirement c^nn2 < 0.7, where wnn2 is 
a neural network output defined for jets containing a secondary vertex. The inputs for 
co>NN2 are similar to the inputs for c^nn, but also the impact parameter significances in the 
x — y and the z plane with respect to the secondary vertex are used. 

From all tracks that pass the previous selection criteria, only the one with the highest 
longitudinal momentum with respect to the jet axis, p™ ax , is kept for each B candidate, 

A reduction of B decay track background in the p™ ax sample is obtained by the require- 
ments cjnn < 0.80 and cunn2 < 0.50 

Fragmentation tracks in the p™ ax sample are removed with the requirement wnn > 0.20. 
Since wnn is designed to achieve optimal separation of b hadron decay tracks from frag- 
mentation tracks using impact parameter information and kinematics, the Monte Carlo 
indicates a fairly flat cunn distribution for signal pions. On the contrary, fragmentation 
tracks peak at zero. 

A momentum of p > 1.0 GeV/c is required for signal pion candidates. In the simulation, 
the momentum distribution of signal pions has a mean value of 2.9 GeV/c with an RMS 
of 1.3 GeV/c before the cut is applied. The B} mass spectrum for single pion transitions 
is not influenced by the momentum requirement. 



n If no secondary vertex is present in the jet, the cut wnn < 0.50 instead of wnn2 < 0.50 is applied. 



15 



6.2 Btt^ mass spectrum 



The signal pion candidate passing the selection cuts described in Section |6J] is combined with 
the corresponding B candidate to form a B} candidate. The invariant mass is calculated using 
Equation |6] as for the B* mass, where the photon is replaced by a pion and the appropriate 
pion mass term is added. The cuts of the signal pion selection have been chosen to obtain an 
acceptable signal-to-background ratio at high signal efficiency. The order of the non-commuting 
selection requirements using c^nn, ^nn2 and p™ ax aims to maximise the difference between the 
shape of the signal and background contributions to the Mb* distribution. 

Due to the intrinsic widths of the B} states and the limited detector resolution, only a 
single peak is seen in the M-q-k spectrum of Figure 0a on top of the combinatorial background. 
According to the simulation, the M-q-k resolution can be described by the sum of a narrow 
Gaussian and an asymmetric Gaussian, both constrained to the same mean value (see Figure |8|). 
The mass resolution depends linearly on M^ n . In the B} signal region around 5.7 GeV/c 2 , the 
standard deviation of the narrow Gaussian is a = 33 MeV/c 2 , and 85% of the resolution function 
entries are contained within 3a. The reliability of the simulated B meson energy and direction 
resolution which dominate the B7r mass resolution is proven by a well simulated shape and peak 
position of the B* signal using the conversion photon sample (see Section [5]). 

The Monte Carlo combinatorial background is checked against data using different test 
samples strongly enhanced in each of the following physics background sources: 1) Fake B} 
candidates from light and charm quark events; 2) Fake B} arising from true b hadrons combined 
with a pion from the weak decay of the b hadron itself; 3) Fake B} formed by combining true b 
hadrons with fragmentation tracks which have not come from a B} resonance. The simulation 
indicates that each test sample is strongly enhanced in the background source under study and 
that the B} signal is suppressed by about a factor of eight compared to the original B} signal 
selection. The Bit mass distributions of the background samples in data are compared with 



the corresponding Monte Carlo mass distributions (see also Section 8.1.2 ). In the case of a 
significant deviation, the simulated background is reweighted to the data. The Monte Carlo 
background distribution so obtained is fitted using a threshold function of the form 



G x -^jx- (m B + m n ) ■ ($ j , (7) 

where $ is the Landau density function^. This background function gives a good empirical 
fit with only four free parameters Cj. The fitted Monte Carlo background is normalised in the 
sideband region 6.10 GeV/c 2 < M^n < 7.10 GeV/c 2 and subtracted from the data distribution. 
The obtained B} signal is shown in Figure 0b. 

The reconstruction efficiency for B} depends on the reconstructed mass M^ n . Monte 
Carlo studies indicate that the efficiency stays constant at high values down to = 
5.7 GeV/c 2 . Below 5.7 GeV/c 2 , the reconstruction efficiency becomes smaller as M^n decreases, 
mainly due to the pp ax requirement. At the B7r mass threshold, the signal efficiency is close 
to zero. The B} distribution of Figure 0b is corrected for efficiency and the resulting signal is 
shown in Figure 0c. 

The mean mass, shape and yield of the observed B} signal is roughly in agreement with 
other measurements The structure of the B^ mass spectrum is too broad to stem from 

a single resonance and leaves room for interpretation. The peak is expected to contain two 
broad and two narrow B} states, and part of the true mass spectrum is shifted to lower mass 

i2$ (A ) = i f c+ioo e Xs+lns d s 

v ' Z7T2 J C — ZOO 
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values by 46 MeV/c 2 due to the omission of the photon in the reconstruction of B} — > B*ir 
decays. The peak may also include a small fraction of B* sJ due to the misidentification of kaons 
as pions. In addition, the peak may contain contributions from B} — > B^nX giving rise to 
satellite peaks in the region 5.4 GeV/c 2 < M-q-k < 5.6 GeV/c 2 , since X is not included in the 
invariant mass calculation. If broad B} states have masses close to the Bn threshold, they 
have an asymmetric signal shape due to phase-space suppression. Radially-excited B mesons 
(2S) in the decay channels B^' -»■ B^tt and B^' -> BW 7T7T may contribute, although the 
production rate of B^' is assumed to be small compared to the B} production rate according 
to ]IU| and P2| . Since there are several ambiguities, e.g. due to B} — > B^irir, B^' decays 



and uncertainties in the combinatorial background, further details of the signal can only be 
obtained by making additional, model-dependent assumptions. 



7 transitions to B* and to B 

In this section, information from the reconstructed Bir^ mass is combined with the weight 
W(B*). The total B} sample is divided into two samples, one enriched and one depleted in 
the decay B} — > B*n(X), by applying a cut on W(B*), as indicated in Figure |a. These two 
B} samples allow a model-independent measurement of the ratio BR(B} — > B*7r(X)), where 
no distinction between decays to B*tt and B*irX is possible. Also a fit with BR(B} — > B*tt) as 
an additional free fit parameter is performed. In this fit, only B} — > B^ir and B} — > B^nn 
decays are considered. The fit result of the branching ratio BR(B} — > B*7r) does not include 
the B*7T7r final state as the sensitivity of the fit is negligible for di-pion transitions of the B}. 



7.1 Model-independent measurement of BR(B} — > B*7r(X)) 

The branching ratio BR(B} — > B*ir(X)) is obtained by counting the number of signal en- 
tries of the B} samples enriched or depleted in the decay B} — > B*tt ± (X). The cut value 
on W(B*) is chosen to minimise the uncertainty of the measurement. The statistical error on 
BR(B} — > B*7r(X)) is minimal if both subsamples are of the same size. Systematic uncertainties 
in the Btc^ background have minimal impact on BR(B} — > B*7r(X)) if the signal-to-background 
ratio is the same for the Btt^ 1 mass distributions of the B*-enriched and the B*-depleted sam- 
ples. The optimal cut on W(B*) is 0.648, fulfilling the minimal systematic error requirement 
and coming as close as possible to the minimum statistical error requirement (see Figure ^). 

The B* enrichment and depletion method can be evaluated by the different selection ef- 
ficiencies for the transitions B} — > B*^ and B} — > B^ in the B*-enriched and B*-depleted 
samples. With the definitions 

• £* : B} — > B*7r efficiency of B*-enriched sample; 

• e* : B} — >• B*7r efficiency of B*-depleted sample; 

• e E : B} — > Bit efficiency of B*-enriched sample; 

• e D : B} — > Bn efficiency of B*-depleted sample; 

• efficiency ratios : e = £ D /s* ; e = e E /e ; e* = e* /e* , 

we calculate from Monte Carlo the efficiency values presented in Table |2|. The numbers reflect 
the cut on W(B*) and thus the quality of the B* enrichment versus the B* depletion. Only the 
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efficiency ratios given in the right column of Table || not the absolute efficiencies, are needed 
for the determination of BR(B} — > B*7r(X)). For the BR(B} — > B*rr(X)) measurement, the 



efficiency 


efficiency ratio 




0.05084 ±0.00023 


e 


1.566 ±0.010 




0.06680 ±0.00026 


e 


0.578 ±0.005 




0.04601 ± 0.00030 


e* 


1.314 ±0.008 




0.07962 ±0.00039 







Table 2: Efficiencies for the reconstruction of B} decaying to B*7r and to B7r. The numbers are 
calculated with respect to the total number of B} passing the B selection in the Monte Carlo. 
Therefore, the numbers reflect the effect of the cut on W(B*) and the charged pion selection. 
Also a factor of 2/3 assuming isospin symmetry to account for decays of B} via neutral pions 
is included in each of the efficiency values. The errors are statistical only. 



invariant Bn^ 1 mass distributions of the B*-enriched and the B*-depleted sample are used. 
Both mass distributions are independent subsamples of the distribution shown in Figure [7|, but 
contain different compositions of B} — > B*tt ± (X) and B} — > Brr^X) decays. Figures |§and [10| 
show the Brr 1 * 1 mass distributions for the B*-enriched and the B*-depleted sample, respectively. 
The Monte Carlo background distributions of both samples are corrected using a procedure 
explained in Sections |6] and |8| and the same fit as for the background of the total B} sample (see 
Section |6|) is performed. The fitted background functions are normalised in the sideband region 
6.10 GeV/c 2 < M Bn < 7.10 GeV/c 2 and subtracted from the corresponding data distributions. 
From the resulting signal peaks, BR(B} — > B*tt(X)) is obtained. Assuming that B} decay to 
B* or B only and with the efficiency ratios defined as above, we derive the following formula: 

BR(B} — » B*t(X)) - BR(B^BMX)) 



BR(B} -> B*rr(X)) + BR(B} -> Btt(X)) 

Ne - e ■ iV D 

60 ' (e -e*)-N E +(l-e-e Q )-N D ' U 

where N E (Nd) denotes the number of B} signal entries of the B*-enriched (B*-depleted) sample. 
In the data, N E = (8782 ± 252 stat ) and N D = (12051 ± 295 stat ) B} candidates are observed in 
the Mb,t signal window (5.3 — 6.1) GeV/c 2 . Using the numbers for the efficiency ratios eo, e 
and e* presented in Table 0, we calculate BR(B} — > B*7r(X)) = 0.85. 

The statistical errors on N E and result in a total error on the branching ratio of ±0.13. 
Besides this error, statistical uncertainties due to the sideband normalisation have been taken 
into account. Since the samples are mutually exclusive, the statistical errors of the sideband 
normalisation of both samples are independent. The contributions of the B*-enriched and the 
B*-depleted sample sideband normalisation to the statistical error on BR(B} — > B*7r(X)) are 
toas an d ±0.15, respectively. Adding all quoted errors in quadrature, the branching ratio of 
orbitally-excited B mesons decaying into B* is measured to be 



BR(B} B*tt(X)) = 0.85 



t-0.26 
-0.27 i 



where the error is statistical only. This branching ratio includes all decays of the type B} — > 
B^rrX, as no cut against additional B} decay products is applied. Consequently, the notation 
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BR(B} — * B*7r(X)) is chosen. Systematic uncertainties, especially of the efficiency ratios and 
the combinatorial Bit background are discussed in Section || 

We further investigate the composition of the B} sample by splitting the sample into 
B} — > B*7r ± (X) and B} — ► B7r ± (X) components. By subtracting from the Btt mass distribution 
of the B*-enriched sample the corresponding distribution of the B*-depleted sample multiplied 
by a scale factor, a B7r ± mass distribution containing B} — > B*tc ± (X) transitions only is ob- 
tained. The scale factor is the ratio of the B} — ■> Btt^ efficiencies of both samples, e = e e /e . In 
a similar way, a mass distribution with B} — ► B*tt ± (X) decays subtracted off is obtained. The 
corresponding efficiency-corrected Bir^ mass distributions for pure B} — > B*tt ± (X) and pure 
B} — > B7r ± (X) transitions are shown in Figure [11]. The number of signal entries in the Bir^ 



mass distributions of Figures 11a and lib depends on the ratio BR(B} — > B*ty(X)) as well as 
on the efficiency ratios defined in Table |2|. 

A significant excess of entries is seen in the pure B} -> B*tt ± (X) distribution at masses 
around 5.7 GeV/c 2 with tails down to 5.5 GeV/c 2 and up to 6.0 GeV/c 2 . The narrow peak in 
the B} -> B*n ± (X) distribution is most likely due to B x (3/2) -> B*tt ± and B* -> B*vr ± decays. 
To obtain the true mass values of the B*7r states, the entries have to be shifted to higher masses 
by 46 MeV/c 2 . 

In the pure B} — > B7r ± (X) mass distribution, a small excess is observed in the region up to 
5.85 GeV/c 2 . This excess can be assigned to the decays B2 — > Btt^ and Bq — > Bn ± . Since the 
statistical significance of the excess in the B} — > B7r ± (X) mass distribution is small, no further 
conclusion is drawn from Figure Ob. 



7.2 Simultaneous fit to the Btt^ mass spectra 

A simultaneous fit is performed to the background subtracted and efficiency corrected Btt mass 
spectra shown in Figures |^c and 10c. Several assumptions are made on the nature of the 
observed signal to reduce the number of free fit parameters: 

• The signal excess stems from B} decays only. Contributions from B* sJ decays, B^' — > 
B^tt and B«' -> B« 7T7T do not exceed a few percent |10|, ^1 and are therefore not 
implemented in the fit but considered as sources of systematic uncertainties. Any other 
excited state eventually contributing to the signal peak is ignored since there is no ex- 
perimental evidence for such states and theoretical predictions give negligible production 
rates. 

• The heavy quark limit itlq — > 00 holds to describe the four B} states. Therefore, according 
to spin-parity conservation, one expects five different mass peaks from single pion tran- 
sitions as listed in the first paragraph of Section |l] and shown in Figure |l[ Furthermore, 
the physical B x states are B^l/2) and Bi(3/2) and thus no mixing occurs. 

• Partners of the same doublet are assumed to have similar properties. The constraints 
on masses, widths and production rates used in the fit are presented in Table [| The 
mass splitting between the narrow states can be calculated using the corresponding mass 
splitting of the D} which has been measured [fj]. The mass splitting between the broad 
states is expected to be of about the same size. Also the order of magnitude of the widths 
of the narrow and broad B} states can be estimated from experimental D} results J 



Only the decay modes listed in Table |3] are taken into account. We explicitly allow the 
decay via two pions to B* and B for all B} states. For each B} state, we set BR(B} — ► 
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state 




fit constraints 


allowed 
decay modes 


production rate 


mass 


width 


Bi(|) 

Ml) 

B* 


°l/2 

1 + 

-■-1/2 

1 + 
L 3/2 

2+ 

^3/2 


f(b - Bg) 
f(b - BS) 
| • f(b - BS) 
| • f(b -> B^ 


free 

M(BS) + 20 MeV/c 2 
free 

M(Bi(|)) + 12 MeV/c 2 


free 

1.25 -r(B5) 

free 

Loo-rCBid)) 


I37T, O 7T7T, r>7T7T 

B*7r, B*7T7r, Btttt 
B*tt, B*tttt, Btttt 
B*tt, Btt, 
B*7T7r, Btttt 



Table 3: Constraints on production rates, masses and widths used in the fit to the total Btt 
mass spectrum. 



B*7T7r) = BR(B} — ► B7T7r) and we assume the same branching ratio BR(B} — > B^tttt) for 
all B} states. For the B*, we set BR(B* -> B*vr) = BR(B* -> Btt). 

The fraction of narrow states f nar row — (f(b — > Bi(3/2))+f(b — > B;!;))/f(b — > B}) is fixed to 
0.6. This number is the average of 1/2, 2/3 and 2/3 corresponding to production rates of 
narrow B} according to state counting, total spin counting and light quark spin counting, 
respectively. To justify this constraint, a fit to the total Btt mass spectrum is performed 
with f narrow as an additional free parameter. 

The mass splitting between B and B* is fixed to the current world average [Q]. 

Each of the five single pion decay modes is represented by a Breit-Wigner function con- 
voluted with the Mbtt dependent resolution function explained in Section |6.2| and shown 
in Figure [8|. To take into account the phase-space suppression at threshold, asymmetric 
Breit-Wigner functions with the threshold factor |35[| 



/threshold (M) 



are used instead of the symmetric Breit-Wigner functions for the broad states. 

• For the di-pion transitions, the signal shape including the detector resolution is taken from 
simulated B} — > B^tttt decays. Simple Gaussians truncated at threshold give a good 
description of the simulated satellite peaks. The mean of the Gaussian depends linearly 
on the mass difference between B} and B ground state. The width of the Gaussian is also 
a function of this mass splitting and depends on the width of the B} state. The functions 
to parameterise the mean and width of the Gaussians are taken from the simulation. No 
attempt is made to implement different signal shapes for decays where the two pions form 
an intermediate resonance or for cascade transitions from high mass B} states via low 
mass B} states to the ground states B* and B 0. 

The different reconstruction efficiencies for B} — > B*^ and B} — > B7r ± decays of the B*- 
enriched and B*-depleted signals are taken from Monte Carlo (see Table 0). As the specific 
peaks that make up the total B} signal have different sizes in the B*-enriched and B*-depleted 

13 Strong decays within the B} multiplet are allowed if the mass splitting within the multiplet or some of the 
widths are larger than the pion mass. 
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mass distributions, different enhancements or depletions are expected for specific regions of 
the B^ mass spectra according to the assumptions on the nature of the B} signal. Thus the 
simultaneous fit provides a non-trivial consistency check of the BR(B} — ► B*tt(X)) result of 
Section |7.1| and of some of the constraints used in the fit. 

The implementation of BR(B} — > B*7r) as a fit parameter causes some complications. Note 
that the BR(B} — ► B*tt(X)) result of Section |7.1| includes transitions of narrow B} and broad 
B} via emission of one pion (and perhaps other decay products). The branching fraction 
to the B* ground state might be different for narrow and broad states and also different for 
B*7r and B*7T7t final states. Whereas the result obtained from Equation ^| is the average of 
the natural composition of the different decay modes, four different parameters have to be 
considered for this fit: BR(B}„ - B*tt), BR(B} broad - B*tt), BR(B}„ - B*tttt) 
and BR(B} broad — > B*7r7r). Monte Carlo studies indicate that the sensitivity of the fit to 
BR(B} narrow — > B*7T7r) and BR(B} broad — > B*tttt) is negligible since shape and position of the 
corresponding peaks in the B71" mass spectra are almost the same for the B*tttt and Bun 
final states. The sensitivity to BR(B} broad — > B*tt) is also smaller than the corresponding 
sensitivity to BR(B} narrow — > B*tt) since a large width but a comparable intra-doublet mass 
splitting of the B} states makes a separation of the decay modes of broad B} to B*tt and Btt 
difficult. To reduce the number of fit parameters and to keep the correlations between the fit 
parameters small, BR(B}„ - B*tttt)/(BR(B}„ - B*tttt) + BR(B}„ - Btttt)) and 
BR(B} broad - B*7T7r)/(BR(B} broad - B*tttt) + BR(B} broad - Btttt)) are fixed to 0.5 in the fit. 
Furthermore, we require BR(B} narrow — > B*tt) = | • BR(B} broad — > B*tt) and the fit parameter 
BR(B} — > B*tt) is the weighted mean of both numbers according to the production rates of 
broad and narrow B}. The factor | is based on the assumption of the same production rates for 
states within the same doublet and BR(B* -> B*tt)/(BR(B* -> B*tt) + BR(B* -> Btt)) = 0.5. 
The fit constraints explained in this paragraph are listed in Table |]. 



fit constraint 


value 


BR(B* ^B*7T7r) 

v -'narrow 


1 


(BR(B* ->B*7r7r)+BR(B* -^Btttt)) 

v v ^narrow ^narrow 


2 


BR(B* ^B*7T7r) 


1 


(BR(B* - > B*7T7r)+BR(B* ^Btttt)) 

''broad J broad 


2 


BR(B* ^B*tt) 

v -'narrow 


3 


BR(B* ^B*tt) 

" broad 


2 



Table 4: Additional fit constraints used in the simultaneous fit. 



With the constraints discussed above (see also Tables |3] and |j) the remaining free parameters 
are: sum of the number of entries of the B*-enriched and the B*-depleted signal, M(Bi(3/2)), 
r(Bi(3/2)), M(B*), T(B*), BR(B} -> B^tttt) and BR(B} -> B*tt). A fit @ is performed to 
the simulated B} signals using the full Monte Carlo statistics. All fit results lie within la of the 
Monte Carlo input value. We fit the Btt ± mass spectra of Figures |9|c and pT0| c simultaneously. 
The least squares fit is performed in the Btt mass region of 5.40 — 6.10 GeV/ c 2 with a bin width 
of 20 MeV/c 2 . The fit results are 

M(B 1 (3/2)) = ( 5.7381^06 )GeV/c 2 
r(Bi(3/2)) = (18^)MeV/c 2 
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BRfB 



.1 

BRfB 



M(B*) 

r(B*) 

BWtttt) 
B*tt) 



(5.839 
( 1291 
0.2451 
0.74 



-0.013 
- 0.014 
-27 

- 23 

f 0.027 

- 0.028 
0.12 
0.10 



GeV/c 2 



MeV/c 2 



(*) 
(*) 



where the errors are of statistical origin only. The fit probability is 65% and the result is 
presented in Figures O and [T3|. The fit results are in agreement with the interpretation of the 
B} signals of Figure [11] given in Section 7.1 . Numbers labelled with (*) should be taken with 
care because of large systematic errors. The robustness of the fit results will be discussed in 
Section |[ 

Note that the statistical error of the BR(B} — > B*rr(X)) measurement presented in Sec- 
tion |7.1| includes the errors arising from the sideband normalisation. The statistical errors of 
the fit results on the other hand do not include the sideband normalisation errors. The latter 
will be discussed in Section ||. The BR(B} — > B*tt) result of the fit does not include decays 
to B* via di-pion emission. The BR(B} — > B^irn) result is corrected by a factor of 0.75 to 
account for double counting of the B^Tr+Tr - final state. In Table the correlations between 
all fit parameters are shown. Systematic uncertainties of the fit results are discussed in detail 
in Section M. 



parameter 


N(B}) 


M(B X ) 


iW 


M(B*) 


r(B*) BR(Bj - 




BR(B} -> B*tt) 


N(B}) 


1.000 


0.078 


0.349 


-0.108 


0.418 


-0.088 


0.114 


Af(Bi(3/2)) 


0.078 


1.000 


0.394 


0.028 


-0.032 


0.067 


0.731 


r(B 1 (3/2)) 


0.349 


0.394 


1.000 


-0.675 


0.355 


-0.764 


0.791 


M(BS) 


-0.108 


0.028 


-0.675 


1.000 


-0.313 


0.741 


-0.380 


r(BS) 


0.418 


-0.032 


0.355 


-0.313 


1.000 


-0.473 


-0.036 


BR(B} bMtttt) 


-0.088 


0.067 


-0.764 


0.741 


-0.473 


1.000 


-0.437 


BR(B} -> B*tt) 


0.114 


0.731 


0.791 


-0.380 


-0.036 


-0.437 


1.000 



Table 5: Table of correlations of all free parameters in the final fit to the data mass spectra of 
the B*-enriched and B*-depleted samples. 



8 Systematic uncertainties 

In the following sections, the determination of the systematic uncertainties is presented sepa- 
rately for the model-independent BR(B} — > B*tt(X)) measurement and for the results obtained 
from the simultaneous fit to the Btt ± mass spectra with different B* content. 

8.1 Systematic error on the BR(Bj — > B*7r(X)) measurement 

For the BR(B} — > B*7r(X)) measurement, the dominant sources of systematic error are un- 
certainties in the efficiency ratios, the modelling of the combinatorial Brr ± background and 
systematic errors on the sideband normalisation of the B*-enriched and B*-depleted samples. 
Each contribution to the total error on BR(B} — > B*tt(X)) is listed in Table |6|. 
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8.1.1 Reconstruction efficiencies 

Monte Carlo simulations are used to calculate the efficiency ratios e, e* and cq. The system- 
atic errors on these ratios are dominated by uncertainties in the photon reconstruction. The 
simulation is checked against data using known properties of B* and 7r°. The latter are formed 
by a pairwise combination of two good conversion candidates or one good conversion and one 
good calorimeter candidate assigned to the same B candidate. 

• We perform direct checks of the photon reconstruction efficiencies: The yields of the Monte 
Carlo B* and 7r° mass peaks are consistent with the results observed in data with the 
simulated B* and tt° production rates being in agreement with earlier measurements @. 
To account for the statistical error of the number of B* and tt° peak entries and for 
possible uncertainties in the simulated production rates, the calculation of efficiency ratios 
is repeated on Monte Carlo with the reconstruction efficiency of conversion (calorimeter) 
photons in the decay B* — > B7 changed by ±10% (lio%). The variation of +15% reflects 
a small discrepancy observed in the simulated and measured B* yields of the calorimeter 
sample. 

• In a B} decay, the helicity angle 6* is the angle between the signal pion momentum 
measured in the B} rest frame and the momentum of the B} in the lab frame. As the 
signal pion selection acceptance depends on cos#*, the B} efficiency is sensitive to the 
shape of cos#*. The distributions of helicity angle for Bq and Bi(l/2) decays are assumed 
to be flat (S-wave transitions) and according to [|37| the Bi(3/2) and the B2 are expected 



to have the same cos 6* distribution for any initial b polarisation: 

(l ;(Bi(3/2),B2 -»• B,B*tt) = - (1 + 3cos 2 #* -6w 3/2 (cos 2 r - -) ) (10) 



rdcosfl* v v ' " £ ' ' 4 V ' v 3 

where W3/2 is the probability that fragmentation leads to a state with the maximum 
helicity value of 3/2 for the light degrees of freedom. The Monte Carlo cos#* distributions 
of Bi(3/2) and B2 have been reweighted to cover the whole range W3/2 = — 1. 

The number of good calorimeter photon candidates iV 7ECAL (Figure ^|b) is not modelled 
well in the simulation. Therefore, the Monte Carlo distribution is reweighted to the 
corresponding data distribution. The reweighting clearly improves the general agreement 
between data and Monte Carlo and has an impact on the efficiency rations e, e* and e$. 
The central value of BR(B} — > B*7r(X)) changes by -0.059 due to the reweighting. To 
quantify the uncertainty in the reweighting procedure we take half of the total change of 
the central value as the error on BR(B} — > B*tt(X)). 

The mass dependent efficiency correction to the B} signal for both the B*-enriched and 
B*-depleted samples produces a deviation in BR(B} — > B*7r(X)) relative to the result 
without mass dependent efficiency correction. To account for any mismodelling in the 
simulated B} masses and the simulated mass dependence of the efficiency, half of this 
deviation is taken as the systematic error. 

The calculated efficiency ratios are uncertain due to limited Monte Carlo statistics. 

We check the calibration of the photon energy measurement by comparing the measured 
and simulated shapes and peak positions of both the B* and the it . There is agreement 
within the statistical errors. This results in small (negligible) errors on the efficiencies of 
the calorimeter (conversion) sample. 
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8.1.2 Background related uncertainties 

Uncertainties in the shape of the simulated background have an impact on the number of signal 
candidates Ne and N-^,. Since the combinatorial backgrounds in both the B*-enriched and the 
B*-depleted samples are affected by systematic shifts in a similar way, the measurement is 
rather robust against possible uncertainties in the Monte Carlo background simulation. For the 
determination of systematic errors, the simulated Bit background is varied using the methods 
described below. For each variation, the Monte Carlo background is normalised and subtracted 
from the data and the number of signal entries Ne and Np are counted. 

• Data test samples are developed in which individual background sources are substantially 
enhanced to study the Bir^ 1 combinatorial background. The background is divided into 
three different classes: tracks combined with mistagged B candidates in light and charm 
quark events (udsc flavour), b hadron decay tracks combined with true b hadrons (b 
hadron decay) and b fragmentation tracks combined with true b hadrons (b fragmenta- 
tion). The selection criteria for each test sample are chosen to cover a large fraction of the 
kinematic region of the signal pion selection. A purity of at least 90% for the background 
source under study and a B} signal fraction smaller than 1.5% is obtained by inverting 
cuts in the original B} selection. For the light and charm quark background, also a sample 



of D* + candidates reconstructed as described in |38[ is used as a qualitative cross check. 
For each test sample, the Bir^ invariant mass distribution observed in data is compared 
with the Monte Carlo distribution normalised to the same number of selected B candi- 
dates. The mass distributions and their bin-by-bin ratios data/Monte Carlo are shown in 
Figure [TJ|. The different ratios are fitted with simple polynomials. The latter are used to 
correct the shape of the original Monte Carlo B^ mass distributions for each background 
source separately. The systematic uncertainty on each background source is given by the 
difference in BR(B} — > B*7r(X)) between using the corrected and the uncorrected shape 
of the B^ mass distribution. 

• The composition of the Bn^ 1 background, as seen in the Monte Carlo after the corrections 
have been applied, is varied for each source. The fraction of each of the three background 
sources described earlier is varied by ±20%. 

• The Peterson fragmentation parameter eb has been varied in the range 0.0028 — 0.0057 
to cover uncertainties in the average fraction of the beam energy carried by the weakly 
decaying b hadron, (xe), and in the shape of the fragmentation function. This variation 
causes a minor change in the Bir^ 1 background shape. The effect on BR(B} — > B*tt(X)) 
is smaller than 0.002. 

• The average charged multiplicity of weakly decaying b hadrons (including K° and A decay 
products) is varied in the range 5.375 — 5.865 in the simulation (see [[39)). The observed 
effect on the Btt^ 1 background shape results in a negligible change of BR(B} — > B*7r(X)). 

8.1.3 Other sources of systematic uncertainties and consistency checks 

The following systematic studies have been performed in addition to the studies described in 
Sections |8.1.1| and |3.1.2| . 



The range of the sideband used for the background normalisation is varied by ±100 MeV/c 2 
on each side for both the B*-enriched and B*-depleted sample. This variation is motivated 
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source 


range 


A(BR(B} -> B*tt(X))) 


B* efficiency (ECAL) variation 

B* efficiency (7 — > e + e~) variation 

cos#* dependency 

reweighting of iV 7ECAL 

Mbtt dependence of B} efficiency 

statistical error on efficiency ratios 

reconstructed B* mass (ECAL) 


[0.90,1-15] 
[0.90,1.10] 

W3/2 e [0, 1] 

« 1% 
±2 MeV/c 2 


±0.036 -0.047 
±0.019 -0.018 
±0.040 -0.034 

±0.030 

±0.018 

±0.018 
±0.007 -0.005 


relative composition of background sources 
B decay tracks background modelling 
b fragmentation tracks background modelling 
udsc tracks background modelling 


±20% 
corr. on/ off 
corr. on/off 
corr. on/off 


±0.027 -0.037 
±0.017 
±0.005 
±0.002 


sideband range variation 
variation of cuts on Vv{B ) 
B* sJ reflections 
B^' reflections 

track parameter resolution variation 


±100 MeV/c 2 
±10% 


±0.076 -0.057 
±0.030 -0.043 
±0.006 -0.026 
±0.000 -0.017 
< 0.010 


total 




±0.12 -0.12 



Table 6: Systematic errors of the BR(B} — > B*7r(X)) measurement. Detailed information for 
each uncertainty is given in the text, as well as a discussion of uncertainties which are negligible 
and thus excluded from this table. 



by the range and shape of the B} signal observed in Figure 0. The quadratic sum of the 
differences observed in the number of signal entries gives the largest error contribution to 
the systematic error of the BR(B} — ► B*ir(X)) measurement. 

The cut on W(B*) has been varied. All cut values producing a ratio of signal to noise 
ratios of the B*-enriched and B*-depleted samples between 0.9 and 1.1 are considered. 
The observed deviations in BR(B} — > B*n(X)) do not exceed to'x^ which is taken as the 
systematic error. 

The amount of B* sJ seen in the B} peaks is less than 4% for a B* sJ production rate 
consistent with ||. The branching ratio BR(B*j — > B*K) is varied from 0.2 to 1.0 in the 
simulation. 

Contributions from radial excitations of B mesons decaying to B^n or B^nn may be 
present in the B} signal [|T0|, [33f] . Monte Carlo studies with simulated B' — > B7r + 7r~ 
and B*' — > B*7r + 7r~ decays^ indicate a contamination of the B} signal around = 
5.6 GeV/c 2 . With the total B^' production rate observed in [p3| , the B} peak of Figure [F] 



does not contain more than 3% B^' transitions. The fraction of B' in the simulated B^' 
sample is varied from 0.3 to 0.7. 

To account for any uncertainties arising from a wrongly simulated tracking resolution, 



the reconstructed track parameters are smeared by ±10% in the Monte Carlo [28 



The cuts of the signal pion selection have been varied. No systematic deviations are 
observed. 



'using M(B') = 5.883 GeV/c 2 and M(B*') = 5.898 GeV/c 2 according to O and in agreement with 
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• The whole analysis is repeated using conversion photons only and calorimeter photons 
only The obtained BR(B} -> B*tt(X)) results and the B} -> B*vr ± (X) and B} -> Btt ± (X) 
mass distributions are in agreement with each other and with the total sample. 

• Varying the cut on £> eve nt so that the b purity changes from 92% to 98% produces no 
systematic deviation in BR(B} — ► B*tt(X)). 

• A neural network has been trained to replace the weight W(B*) obtained from the Monte 
Carlo purity parameterisation. The neural network output is strongly correlated with the 
weight W(B*) and does not improve the B*/B separation. 

All systematic errors considered for BR(B} — ► B*rr(X)) are added up in quadrature, result- 
ing in a total error on BR(B} B*tt(X)) of ±0.12 (see Table §). 

8.2 Systematic uncertainties on the results of the simultaneous fit 
to the B*-enriched and B*-depleted mass spectra 

All sources of systematic error are varied in turn within the estimated uncertainty range and 
the fit presented in Section [7^ is repeated. If not stated otherwise, the observed deviation 
with respect to the original fit is taken as a systematic error for each fit parameter. The total 
systematic error on each fit parameter is the quadratic sum of all individual error contributions. 
This procedure takes into account the correlations between the fit parameters. All systematic 
uncertainties are listed in Table 0. 

8.2.1 Variation of fit constraints 

The fit constraints are varied according to the present knowledge of D} properties @ and 
theoretical considerations P,fi|,|TT|-|n|: 

• BR(B} broad — > B*7r) was fixed to | • BR(B} narrow — > B*7r). This constraint is changed to 
BR(B} broad — > B*7r) = BR(B} narrow — > B*ir) and the corresponding deviations of the fit 
parameters are taken as systematic errors. 

• The ratio of the widths of partners of the same doublet was fixed to r(B2)/r(Bi(3/2)) = 
1.25 and r(B 1 (l/2))/r(B^) = 1.00. We allow a variation of the ratio of widths of 1.0-1.4 
for the narrow states and 0.7-1.4 for the broad states. 

• The relative production rates f(b -> B*) : f(b -> Bi(l/2)) : f(b -> Bi(3/2)) : f(b -> B*) 
were fixed to 2 : 2 : 3 : 3. We vary the production rate ratio f(b — > Bi(3/2))/(f(b — > 
Bi(3/2))+f(b -> B*)) in the range 0.375-0.600, f(b -> B*)/(f(b -> B*)+f(b -> Bi(l/2)) in 
the range 0.250-0.600 and f na rrow in the range 0.5-0.75. These variations cover the different 
production rate estimates of state counting (1:1:1:1) to total spin counting (1:3:3:5). 

• The constraint BR(B2 — > B*7r) = BR(B2 — * Bn) was used. We allow a variation of 
BR(B* -> B*tt)/(BR(B* -> B*tt) + BR(B* -> Btt)) = 0.3 - 0.7. 

• Whereas the mass splitting between the j q = 1/2 and j q = 3/2 doublets was free in the 
fit, the mass splittings within the doublets were fixed. We allow a variation of M(B|) — 
M(Bi(3/2)) = (5 - 20) MeV/c 2 and M(Bi(l/2)) - M(B*) = (0 - 50) MeV/c 2 . 
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• The ratio BR(B} -> B*7T7t)/(BR(B} -> B*7T7t) + BR(B} -> Btttt)), which was fixed to 0.5, 
is varied in the range 0.3-0.7. No significant deviation in the fit parameters is observed, 
since the position and shape of the peaks corresponding to the two final states are very 
similar. 



8.2.2 Reconstruction efficiencies 

The B} reconstruction efficiency is a function of the reconstructed B7r mass and the angular 
distribution of the 7T. The efficiency is taken from the Monte Carlo and possible problems with 
its simulation are taken into account: 



The signal pion selection cuts have been varied. Whereas the Monte Carlo simulation 
describes the reconstruction efficiency well at high Bit masses, this statement can not 
be proven for low Bn masses. Therefore, we perform the fit to the B} signal without 
an acceptance correction of the mass spectrum as shown in Figure 0b. Although the fit 
result has a low fit probability, we assign half of the total deviations of the fit parameters 
observed with respect to the original fit as systematic errors. This results in a large 
systematic error on BR(B} — > B^nn). 

Systematic uncertainties arising from the photon reconstruction as well as from the un- 



known cos 9* distributions are determined using the methods presented in 8.1.1 



8.2.3 Background related uncertainties 

• The uncertainties due to the limited data statistics in the upper sideband regions of 
the B*-enriched and B*-depleted samples produce quite large error contributions on 
BR(B} — ► B*7r). For the fit results, these error contributions are treated as system- 
atic errors and have to be compared with the corresponding errors obtained for the 
model-independent BR(B} — ► B*tt(X)) measurement. While the sideband range used 
for the background normalisation is the same for both measurements, the error of the 
BR(B} — > B*7r) fit result is smaller because the fit makes explicit use of the shape and of 
the composition of the B} mass distributions. 

• The Peterson fragmentation parameter and the charged particle multiplicity of weakly 
decaying b hadrons are varied in the simulation as described in Section |3.1.2| . The mul- 
tiplicity variation has no effect on the fit results. 

• To evaluate the uncertainty due to the function used to fit the Monte Carlo background, 
we directly subtract the corrected Monte Carlo background histograms of the Bti mass 
distributions from the corresponding data histograms. 

• The range of the sideband used for the background normalisation is varied by ±100 MeV/c 2 
on each side and the observed deviations are added in quadrature. 



For the determination of the shape and the composition of the Bti background, the same 
methods as presented in Section |S.1.2| are used. 
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source 


ABR 

(Bj - B*tt) 


AM(B 1 (|)) 


Ar(B 1 (D) 


AM (Bo) 


Ar(B^) 


<I75- AB * 

(B* -+ bMtttt) 


BR(B: /broad -B*7r) 

r(B5)/r(Bi(3/2)) 
r(Bi(i/2))/r(B5) 
f 

n&rrow 

prod, rate Bi(3/2) vs. B* 
prod, rate Bg vs. Bi(l/2) 
BR(B^ -» B*7r/B7r) 
M{B* 2 ) - M(Bi(3/2)) 
M(Bi(l/2)) - M(BS) 


+0.059 

+0.014 
-0.009 

+0.001 

+0.049 

— U.Uzo 

±0.012 

+0.004 
-0.003 

+0.013 
-0.008 

+0.003 

— u.yjyjz 

±0.001 


-0.0015 

+0.0008 
-0.0013 

+0.0004 

+0.0001 
— U.UUUz 

+0.0015 
-0.0021 

+0.0001 
-0.0002 

+0.0012 
-0.0004 

±0.0024 
±0.0001 


-4.8 

+4.1 
-1.9 

+0.5 

+7.9 
— o.y 

±1.2 
±0.0 

+2.4 
-1.2 

+2.4 
— o.Z 

+0.0 
-0.2 


+0.0097 

+0.0001 
-0.0029 

-0.0030 

+0.0144 
n ni on 

— u.uiyy 

+0.0000 
-0.0022 

+0.0060 
-0.0054 

+0.0043 
-0.0075 

+0.0017 

— U.UU1U 

+0.0119 
-0.0200 


+1 

+2 
-0 

-20 
±12 

±1 

+i 

-2 

+3 
-2 

±2 

+5 
-6 


+0.011 

+0.007 
-0.016 

-0.004 

+0.001 

— U.UUo 

+0.003 
-0.006 

+0.003 
-0.002 

+0.012 
-0.017 

+0.003 

— \J.\jyjz 

+0.007 
-0.012 


efficiency /(M(Btt)) 
reweighting of JV 7ECAL 
B* efficiency (ECAL) 
efficiency /(cos#*) 


±0.039 
±0.036 

+0.013 
-0.023 
+0.002 
-0.001 


±0.0018 
±0.0021 

+0.0007 
-0.0005 

+0.0001 
-0.0002 


±4.3 
±7.4 

+1.3 
-2.3 
+4.1 

-5.0 


±0.0050 
±0.0025 

+0.0004 
-0.0016 
+0.0077 
-0.0103 


±1 

±2 
+1 

+6 
-7 


±0.109 
±0.021 
±0.003 

+0.002 
-0.004 


sideband norm. B*-enr. 
sideband norm. B*-dep. 
Peterson fragmentation 
background fit function 
sideband range variation 
B decay bg fraction 
B decay bg modelling 
b fragm. bg fraction 
b fragm. bg modelling 
udsc bg fraction 


+0.125 
-0.071 

+0.116 
-0.067 
+0.034 
-0.049 

±0.023 

+0.032 
—0.030 
+0.037 
-0.044 

±0.016 

+0.030 
-0.028 

±0.025 
±0.023 


+0.0041 
-0.0032 

+0.0031 
-0.0025 

+0.0010 
-0.0008 

±0.0025 

+0.0012 
—0.0011 
+0.0001 
-0.0008 

±0.0012 

+0.0003 
-0.0004 

±0.0002 

+0.0001 
-0.0000 


+15.0 
-9.3 

+9.6 
-4.5 

+12.3 
-6.4 

±0.6 

+3.9 
—5.2 

+4.7 
-5.0 

±6.9 

+5.7 
-4.4 

±4.7 

±0.5 


+0.0013 
-0.0020 

+0.0075 
-0.0081 
+0.0044 
-0.0047 

±0.0055 

+0.0037 

— U.0064 

+0.0089 
-0.0082 

±0.0111 

+0.0060 
-0.0084 

±0.0104 
±0.0007 


+3 
-6 

+4 
-9 
+35 
-27 

±6 

+5 
— 13 
+4 
-6 

±25 

+10 

-8 

±33 
±0 


+0.017 
-0.010 

+0.014 
-0.022 

+0.046 
-0.065 

±0.005 

+0.005 
—0.007 

±0.054 
±0.002 

+0.042 
-0.061 

±0.007 
±0.002 


variation of bin width 
a y ' contamination 
B*j contamination 
variation of fit range 


+0.009 
-0.005 

+0.026 
-0.028 

±0.021 

±0.002 


+0.0010 
-0.0008 

+0.0003 
-0.0002 

+0.0002 
-0.0004 

+0.0001 
-0.0000 


+5.0 
-3.5 

+3.5 
-3.7 

±2.7 

+2.0 
-1.7 


+0.0044 
-0.0076 

+0.0055 
-0.0044 

+0.0042 
-0.0043 

+0.0009 
-0.0006 


+5 
-16 
+4 

-5 

±4 

+n 

-10 


+0.017 
-0.033 

±0.003 

+0.004 
-0.005 

+0.010 
-0.006 


total 


+0.21 
-0.15 


±0.007 


+29 

-23 


+0.034 
-0.042 


±63 


+0.143 
-0.161 


central value of fit result 


0.74 


5.738 


18 


5.839 


129 


0.327 



Table 7: Systematic errors on the fit parameters of the simultaneous fit. The numbers for 
ABR(B} — > B^*^7T7r) have to be multiplied by 0.75 to account for double counting. The total 
systematic error of each fit parameter is the quadratic sum of the individual errors. Sources of 
systematic uncertainty with an error of smaller than 10% of the total systematic error for each 
fit parameter are excluded from the table. 
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8.2.4 Contamination of the Bj signal, other uncertainties and cross checks 

• Contributions from higher orbital or radial excitations of B mesons decaying to B^n 
or B^tttt may be contained in the signal peaks of the B*-enriched and the B*-depleted 
samples. The procedure described in Section |8.1.1| is repeated. 



Reflections from B* sJ decays influence the fit results since a small fraction of kaons are 
misidentified as pions. The B* sJ production rate is varied by ±50% to cover the experi- 
mental error ||. 

The range of the fit region is changed by ±80 MeV/ c 2 on both sides. 

The whole analysis is performed using bin widths of 16 and 25 MeV/c 2 instead of 
20 MeV/c 2 . 

A fit to total Bvr ± mass spectrum produces results for M(B 1 (3/2)), r(Bi(3/2)), M(B* ), 
r(B*) and BR(B} -> B^tttt) in good agreement with the fit results obtained from the 
mass spectra of the B*-enriched and B*-depleted samples. Furthermore, the fit to the 
total B7r mass spectrum is repeated with f narr0 w as an additional fit parameter. The result 
fnarrow = 0.76 "t oil stat * s i n agreement with the fixed value of 0.6 used in the original fit. 

The number of B} signal entries observed in Figure [7| and the Monte Carlo reconstruction 
efficiencies produce a total B} production rate consistent with our measurement presented 
in @- 



8.2.5 Robustness of fit results 

Extensive systematic studies have been performed to test the stability of the fit results. Since 
the systematic error on BR(B} — > B^'tttt) is large, the existence of B} — > B^ittt decays can 
be questioned and the fit is repeated with BR(B} — > B^nir) = 0. Furthermore, the functional 
form of the broad B} states, which is not precisely known at Btt threshold, has been varied 
assuming different theoretical approaches. The fit probability for any fit with BR( B} — > B^nn) 
fixed to zero is always below 4% (to be compared with a fit probability of 65% of the original 
fit). Depending on the functional form of the broad B} states, the mass of the broad B} states 
lies below or above the Bi(3/2) mass for BR(B} — > B^nn) = 0. Because of this ambiguity, 
the widths and the mass of the Bq as well as BR(B} — ► B(*^7T7r) are not quoted as robust fit 
results. The mass and the width of the Bi(3/2) as well as BR(B} — > B*n) stay stable for all 
fits. The described ambiguity is also observed in the corresponding fit to the total Bit mass 
spectrum. In comparison to other measurements, especially JT0], this is a rather conservative 
treatment of the fit results. Experiments with comparable mass resolution are expected to have 
similar problems in determining B} properties from a fit to a single Btt mass peak. 



9 Summary and conclusion 

We have analysed orbitally-excited mesons by forming combinations of inclusively reconstructed 
B mesons and charged pions. A new way to determine the combinatorial Btc background using 
data test samples while maintaining high statistics is presented. A high statistics tag of the 
decay B* — > B7 is used to obtain B} samples enriched or depleted in their B* content. We 
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present the first measurement of the branching ratio of orbitally-excited B mesons decaying 
into B*. The result is 

BR(B} -> B*tt(X)) = OMtoir ± °- 12 , 

where the first error is statistical and the second systematic. The measurement does not 
depend on the shape of B7r mass distributions or on any specific model. It is in agreement with 
theoretical predictions and the measured B* and B} production rates at LEP. 

Making further use of the B* information, a simultaneous fit to the B} mass spectra of 
samples enriched or depleted in their B* content is performed. In this fit, the masses, widths 
and production rates of the B 1 (l/2) and B2 are constrained by the corresponding properties of 
their doublet partners Bq and Bi(3/2), respectively. The fit yields 

M(Bi(3/2)) = ( 5.738i°;°^± 0.007 ) GeV/c 2 
T(B 1 (3/2)) = ( 18±S±g )MeV/c 2 
BR(B^B*tt) = 0.741^1^. 

The first error indicates the statistical and the second error the systematic uncertainty. The 
fit favours a contribution of B} — > B^nn decays to the B} signal and a mass of the broad 
B} states about 100 MeV/c 2 above the narrow B} states. Systematic uncertainties in the re- 
construction efficiency and the combinatorial background at low Btt masses together with the 
lack of knowledge of the exact functional form of the broad B} states at Btc threshold do not 
allow an unambiguous determination of the width and mass of the Bg (or Bi(l/2)). The fit 
results are in agreement with predictions from several HQET models. The M(Bj.(3/2)) result 
agrees well with a measurement of M{B^) [§]. On the other hand, a recent B} analysis |10| pre- 
sented masses that disagree with the results of this analysis and with ||. The measured value of 
BR(B} -> BMtttt) = 0.25±0.02l is consistent with the range 0.1-0.2 predicted by theory £| 
and in agreement with an experimental result obtained from the reconstruction of the B(*)7r7r 
final state @. The results for BR(B} ->■ B*n(X)), BR(B} -> B*vr) and BR(B} -> BWtttt) 



are in good agreement with each other on the assumption that the decay channels with X 7^ n 
contributing to BR(B} — > B*7r(X)) are small. 
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B* Jq= 1/2 Jq= 3/2 

A I5 K 2 




electromagnetic: AAAAAAAA 



Figure 1: The four B} states and their dominant decays to the ground state doublet (B, B*). 
Strong decays via single pion emission are indicated as solid (D-wave) and dashed (S-wave) 
lines. The B* decays radiatively because of the small B*-B mass splitting. 
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Figure 2: a) The track weight uj vtx for all tracks in hemispheres with a good secondary vertex. 
The peaks near and 1 correspond to tracks created by b fragmentation and b hadron decay 
tracks, respectively. The peak near 0.5 is produced by tracks which are not unambiguously 
assigned to the primary or the secondary vertex, as in the case of tracks matching both the 
primary and secondary vertex or matching no vertex at all. b) The track weight cjnn for tracks 
of all hemispheres (with or without a good secondary vertex). The separation power of cg>nn 
is superior to the separation power of u vtx . c) The combined track weight cu tr calculated from 
u vtx and cunn for tracks of all hemispheres. Note that cu tT is shown on a logarithmic scale. 
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Figure 3: a) The number of good conversion photon candidates per B candidate observed in 
data and Monte Carlo, b) The number of good calorimeter photon candidates per B candidate 
observed in data and Monte Carlo. For the analysis, the Monte Carlo distribution of the latter 
is reweighted to the data distribution. 
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Figure 4: a) The AM = Mb 7 — Mb mass distribution of the conversion photon sample. The 
background is estimated from Monte Carlo simulation and normalised to the data distribution 
in the sideband region 0.09 GeV/c 2 < AM < 0.20 GeV/c 2 . b) The corresponding background 
subtracted signal. The fit function used for the signal is described in the text. 
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Figure 5: a) The AM = Mb 7 — Mb mass distribution of photons reconstructed in the elec- 
tromagnetic calorimeter. The background is estimated from the Monte Carlo simulation and 
normalised to the data distribution in the sideband region 0.10 GeV/c 2 < AM < 0.20 GeV/c 2 . 
Although the resolution is poor compared to the conversion photon sample, an excess of entries 
in the data distribution around 46 MeV/c 2 is clearly visible, b) The corresponding background 
subtracted signal. The fit function is described in the text. 
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Figure 6: a) The W(B*) distribution for data with the corresponding Monte Carlo histograms 
indicating the number of B candidates with a B* parent and no B* parent. The dotted line gives 
the boundary between the B*-enriched and B*-depleted samples, b) The ratio of the efficiency 
to reconstruct a B meson with a B* parent over the efficiency to reconstruct a B meson without 
a B* parent versus the weight W(B*) calculated from simulated data. 
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Figure 7: a) The B7r ± mass distribution for data. The shaded histogram indicates a fit to the 



corrected Monte Carlo background using a reweighting method described in Section |8.1.2| . The 
function used for the fit is described in the text, b) The B} signal after subtraction of the 
simulated background, c) The efficiency-corrected B} signal. The observed structure of the 
B} signal suggests a superposition of several different states. The mass dependent efficiency 
correction has a strong impact on the signal shape at low Bit mass values. 
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Figure 8: a) Monte Carlo M-q^ resolution of B} decaying to B^it^ 1 in the mass region 
5.65 GeV/c 2 < M(Btt) < 5.70 GeV/c 2 . The fit function is the sum of two Gaussians both 
constrained to the same mean value. Sigmal is the standard deviation of the narrow Gaussian 
and sigma2 (sigma3) corresponds to the left (right) standard deviation of the asymmetric broad 
Gaussian, b) The linear dependence of the width of the resolution function on M^ n is shown 
for each sigma in the B} signal region. 
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Figure 9: a) The ~Bn ± mass distribution of the sample enriched in the decay B} — > B*7r ± (X) in 
data. The shaded histogram indicates a fit to the corrected Monte Carlo background using a 
reweighting method described in Section |8.1.2| . b) The signal distribution after subtraction of 
the simulated background, c) The efficiency-corrected signal. 
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Figure 10: a) The B7r ± mass distribution of the sample depleted in the decay B} — ■> B*7r ± (X) 
in data. The shaded histogram indicates a fit to the corrected Monte Carlo background using 
a reweighting method described in Section g.l.2| . b) The signal distribution after subtraction 
of the simulated background, c) The efficiency-corrected signal. 
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Figure 11: a) The efficiency-corrected Btt 1 * 1 mass distribution of B} — ► B*7r ± (X) transitions seen 
in data. A clear peak is visible at 5.7 GeV/c 2 . The structure is unlikely to stem from a single 
state, b) The efficiency-corrected Bn^ 1 mass distribution of B} — > B7r ± (A) transitions seen in 
data. A 2.2a excess is observed around 5.8 GeV/c 2 . 
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Figure 12: Simultaneous fit to the B7r mass distributions of the B} samples enriched or depleted 



inB} 



B*7r (it) decays (see also Figure [13]). The fit results of the decays of the broad and 



narrow B} for transitions via one and two pions are presented separately. 
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Figure 13: Simultaneous fit to the Bit mass distribution of the B} samples enriched or depleted 
in B} — > B*7r ± (7r) decays. The fit results of the decays of the broad and narrow B} for transitions 



via one and two pions are presented separately. In comparison to Figure [12] the fraction of B} 
decays to B* is reduced with respect to B} decays to B (e.g. compare the light (B} — > B*tyty) 
and the dark (B} — > Btttt) solid lines in Figures [12] and | 
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Figure 14: The Btt^ 1 mass distributions of data and Monte Carlo for each of the three test 
samples (left side) and the corresponding bin-by-bin ratio of the mass distributions (right side). 
a)+b) Light and charm quark sample, c)+d) b hadron decay sample, e)+f) b fragmentation 
sample. 



